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ABSTRACT 

Context. Magnetic heiicity has the remarkable property of being a conserved quantity of ideal magnetohydrodynamics (MHD). 
Therefore, it could be used as an effective tracer of the magnetic field evolution of magnetised plasmas. 

Aims. Theoretical estimations indicate that magnetic heiicity is also essentially conserved with non-ideal MHD processes, e.g. mag¬ 
netic reconnection. This conjecture has however been barely tested, either experimentally or numerically. Thanks to recent advances 
in magnetic heiicity estimation methods, it is now possible to test numerically its dissipation level in general three-dimensional 
datasets. 

Methods. We first revisit the general formulation of the temporal variation of relative magnetic heiicity on a fully bounded volume 
when no hypothesis on the gauge is made. We introduce a method to precisely estimate its dissipation independently of the type of 
non-ideal MHD processes occurring. In a solar-like eruptive event simulation, using different gauges, we compare its estimation in a 
finite volume with its time-integrated flux through the boundaries, hence testing the conservation and dissipation of heiicity. 

Results. We provide an upper bound of the real dissipation of magnetic heiicity: It is quasi-null during the quasi-ideal MHD phase. 
Even when magnetic reconnection is acting the relative dissipation of magnetic heiicity is also very small (< 2.2%), in particular 
compared to the relative dissipation of magnetic energy (> 30 times larger). We finally illustrate how the helicity-flux terms involving 
velocity components are gauge dependent, hence limiting their physical meaning. 

Conclusions. Our study paves the way for more extended and diverse tests of the magnetic heiicity conservation properties. Our 
study confirms the central role that heiicity can play in the study of MHD plasmas. For instance, the conservation of heiicity can be 
used to track the evolution of solar magnetic fields, from its formation in the solar interior until their detection as magnetic cloud in 
the interplanetary space. 
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1. Introduction 


In physics, conservation principle have driven the understand¬ 
ing of observed phenomena. Exact and even approximately con¬ 
served quantities have allowed to better describe and predict the 
behaviour of physical systems. Conservation laws state that, for 
an isolated system, a particular measurable scalar quantity does 
not change as the system evolves. A corollary is that for a non 
isolated system, a conserved scalar quantity only evolves thanks 
to the flux of that quantity through the studied system bound¬ 
aries. Given a physical paradigm, a physical quantity may not 
be conserved if source or dissipation terms exist. 


In the magnetohydrodynamics (MHD) framework, a quan¬ 
tity has received increa sing attention for its conservation prop¬ 
erty: magnetic heiicity (lElsasseijll956l) . Magnetic heiicity quan¬ 
titatively describes the geometrical degree of twist, s hear, or 
more generally, knottedness of magnetic field lines dMoffatd 
Il969l) . In ideal MHD, where magnetic field can be described 
as the collection of individual magnet ic field lines, m agnetic he¬ 
iicity is a strictly conserved quantity (IWoltieii|1958l) as no dissi¬ 
pation, nor creation, of heiicity is permitted since magnetic field 
line cannot reconnect. 


In his seminal work. I^vlorl d 1 974lf conjectured that even in 
non-ideal MHD, the dissipation of magnetic heiicity should be 
relatively weak, and hence, that magnetic heiicity should be con¬ 
served. This could be theoretically explained by the inverse- 
cascade property of magnetic heiicity: in turbulent medium, 
heiicity unlike magnetic energy, tends to cascade towards the 
larger spatial scales , thus avoiding dissip ation at smaller scales 
dFrisch et al.l Il975t iPouauet et all Il97^ . This cascade has 
been observed in numerical simulations dAlexakis et al. 200d 
Mininnil l2007l) as well as in labora tory experimen ts ( .li et al.l 


1995I) . In the case of resistive MHD, iBerger d 19841) derived an 


upper limit on the amount of magnetic heiicity that could be dis¬ 
sipated through constant resistivity. He showed that the typical 
heiicity dissipation time in the solar corona was far exceeding 
the one for magnetic energy dissipation. 


From the Taylor’s conjecture on heiicity conservation have 
been derived multiple important consequences for the dynam¬ 
ics of plasma systems. Based on heiicity conservation. iTavlcM 
dl974l) predicted that relaxing MHD systems should reach 
a linear force free state. This prediction, which was veri¬ 
fied to different degrees, has allowed to understand t he dy- 
namics of plasma in several laboratory experiments dTavloil 
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Il986t IPrageJ[TQQOt Ijil ll999t lYamadal[l999h . The concept of 
Taylor relaxation has also driven theoretic al models of solar 
coronal heating (iHevvaerts & Prie^ Il984ll . The importance 
of magnetic helicity cons ervation has been further raised for 
magnetic field dynamos ([Brandenburg & Subramani^ 120051: 

[Blackman & Hubbardll2014l) . Magnetic helicity is also strongly 

impacting the energy budget during reconnection events 
( Linton et al.ll200l1; iLinton & Antioch^l2002t iDel Sordo et alJ 
l^ldl) . and models of eruption based on magnetic helicity an¬ 
nihilation have been developed (iKusano et al.ll2004ll . The con¬ 
servation of magnetic helicity has been suggested as the core 
reason of the existence of Coronal Mass Ejections (CMEs), the 


latter bei ng the mean for the Su n to expel its excess magnetic 
helicity (lRusti[r9^ lLowlll996l) . Because of this hypothesis, 
important efforts to estimate the magnetic hel icity in the solar 
coronal have been carried over the last decades (lDemoulinl2007l: 
iDemoulin & Pari^l2009h . 

Despite its potential importance, it is surprising to note that 
Taylor’s conjecture on helicity conservation has been barely 
tested. Several experimental and numerical works have focused 
on testing Tay lor’s prediction that relaxed systems should be lin¬ 
ear force free (Tavloi 1986 : [Pra^ll999l :lJil ll999t|Yamadall9^ 


lAntiochos & DeV” I199S ). Howeyer, the fact that the system 

may not reach a complete relaxed linear force free state, does 
not mean that helicity is not well conseryed. It simply means 
that the dynamics of the relaxation does not allow the full redis¬ 
tribution of helicity towards the largest ayailable scales. Looking 
at the long term eyolution of turbulent MHD systems, magnetic 
helicity has been found to decay slower than magnetic energy 
( Biskamp & MulleJl 19991 : IChristensson et al.1l2005t ICandelare^ 


I2OI2I) . These studies are howeyer based on the estimation of the 

helicity density which is not generally gauge inyariant. The peri¬ 
odic boundaries, the long term and large spatial scales inyolyed 
do not allow an estimation of helicity dissipation when a "singu¬ 
lar" non-ideal eyent is occurring. 

Direct tests on magnetic helicity conseryation and dissipa¬ 
tion haye actually been limited because of the inherent difficul¬ 
ties to measure that yery quantity. Laboratory experiments and 
obseryations requires the measurement of the full 3D distribu¬ 
tion of the magnetic held. Most laboratory experiments there¬ 
fore make assumptions on the symmetries of the system to limit 
the sampling o f the data, hence limiting measurement precision. 
iJi et al.l([l995h haye measured a 1.3 -5.1% decay of helicity (rel- 
atiyely to a 4 - 10% decrease of energy) duri ng a sawtooth relax¬ 
ation i n a reyersed-held pinch experiment. iHeidbrink & Dan3 
( I 2 OOOI) haye estimated that helicity was conseryed within 1% 
during a sawtooth crash. Other experiments haye tried to mea¬ 
sure helicity conseryation by comparing the helicity in the sys¬ 
tem with its theo retically injected amount (iBarnes et al.lll98OT 
iGray et akl I 2 OIOI) . While they found results agreeing within 
10 - 20%, the experimental c onditions limit the pr ecision on the 
measure of injected helicity (IStallard et al.ll200^ . In numeri¬ 
cal experim ents of flux tubes reconnect ion, using triple periodic 
boundaries. iLinton & AntiochosI (l2005l) haye found that the loss 
of helicity was ranging between 6% and 53% depending on the 
Lundquist number, hence primarily due to diffusion rather than 
magnetic reconnection. 

Eor non-periodic systems where the magnetic field is thread¬ 
ing the domain boundary (as in most natural cases), the yery def¬ 
inition of magnetic helicity is not gauge inyariant (cf. Sect. EB, 
and a modihed dehnition of magnetic helicity, the relatiye 
magnetic helicity, had to be introduced (iBerger & Eieldl 0*984 
see Sect. lO l. Only yery recently, practical methods that 
allow to generally computed relatiye magnetic helicity haye 


been published ( Rudenko^Myshyakoy 201 it iThalmann et al.l 
I 2 OI lb l^lori et al.ll2012t Yang et al.ll2013[) . These methods are 
now opening the possibilities to carefully study helicity of 3D 
datasets su c h as t he o nes frequently used fo r natural plasmas. 
IZhao et ^ (1201 5h and iKnizhnik et al.l (12015l) haye obseryed a 
good h elicity conseryation but haye not quantihed it. lYang et ^ 
(1201 3h haye measured the helicity conseryation in a numerical 
simulation by comparing the relatiye helicity flux with the yaria- 
tion of helicity in the domain. They found that during the quasi¬ 
ideal eyolution of the system the helicity was conseryed within 
3%, with higher dissipation yalues when non-ideal MHD effects 
become important. Howeyer, in preyious studies, the relatiye 
helicity flux dehnition used may not be fully consistent with 
the choice of gauge used to compute the yolume helicity (see 
Sect. lZ2l) . 

In the present manuscript, we will further push the test on 
the conseryation of magnetic helicity. We will deriye a general¬ 
ized analytical formula for the Hux of relatiye magnetic helicity 
(Sect. I 2 . 3 I 1 without taking any assumption on the gauges of the 
studied and the reference helds. We will discuss whether rela¬ 
tiye magnetic helicity can be considered as a conseryed quan¬ 
tity in the classical sense, that is, whether its yariations can be 
described solely as a hux through the boundary. The general 
method that we will employ (see Sect. [3]) is based on the com¬ 
parison of the eyolution of the relatiye magnetic helicity with its 
Hux through the boundaries, and will be applied to a numerical 
simulation of solar actiye-like eyents. In Sections |4] & |5j using 
different gauges we will constraint the leyel of conseryation of 
relatiye magnetic helicity, and study the amount of dissipated 
magnetic helicity. This is further extended in Appendix lAl with 
another selection of gauges. We will hnally conclude in Sect.|6] 


2. Magnetic helicity and its time variation 

2.1. Magnetic helicity 
Magnetic helicity is dehned as 

f A Bd^, (1) 

J'V 


where B is the magnetic held studied oyer a hxed yolume 'V. 
W is here a fully closecQ yolume bounded by the surface d'V . 
The yector potential, A of B classically yerihes: V x A = B as 
V ■ B = 0 (only ap proximately true in numerical computations, 
I Valori et ^ 120 1 3l) . The magnetic held B is gauge inyariant, 
meaning that it is unchanged by transformations A ^ A -(- Vt/r, 
where if/ is any sufficiently regular, scalar function of space and 
time. Since A is not uniquely dehned in general, magnetic he¬ 
licity requires additional constraints in order to be well-dehned. 
In particular, is a gauge-inyariant quantity proyided that 'V is 
a magnetic yolume, i.e. that the magnetic held is tangent at any 
point of the surface boundary (9'T' of 'T': (B ■ dS) \s^ - 0, at any 
time. 

Assuming that the yolume 'V is hxed, with cl'T' a Hux surface 
ensuring gauge inyariance, the temporal yariation of magnetic 
helicity is deriyed by direct differentiation of Eq. ([TJ: 


df 


f “.Bd-V+ r A.^dT, 

jn/ ot jn/ dt 


( 2 ) 


' We do not consider the particular cases of inhnite volumes which 
require specihc hypothesis on the behaviour of the quantities at large 
distance. 
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where each of the two integrals on the right hand side are well 
defined, since they are independently gauge invariant. Given that 
V ■ (A X dhldt) - dAldt ■ V x A - A • 5(V x A)/dt and using the 
Gauss divergence theorem, one obtains: 


the normal component on the surface B ■ dS Is^y = of the 

studied domain, at any instant there is a unique potential field 
which satisfies the following condition on the whole boundary 
of the volume considered: 



Here, dS is the elementary surface vector, directed outside of the 
domain 'V. Using the Faraday’s law of induction, one derives: 


df/m 

df 




(3) 


Using again the Gauss divergence theorem, one finds that the 
temporal variation of magnetic helicity is composed of three 
independently gauge-invariant terms, Le., a volume dissipative 
term and two helicity flux terms on the surface of "y, such that: 


d/frr 

df 




df 


diss 


Fm,B 



(4) 

(5) 

(6) 

(7) 


In ideal MHD, where E = -v x B, the volume term is null. 
For an isolated system, magnetic helicity is thus conserved in the 
classical sense since its variations are null (cf. Sect. [T]). Varia¬ 
tions of i/m can only originate from advection of helicity through 
the boundaries of ’V. The di/m/dfldiss term corresponds to the 
diss ipation of mag netic helicity of the studied magnetic field in 
"V. iTavlod (119741) conjectured that this term is relatively small 
even when non-ideal MHD processes are developing, e.g., when 
magnetic reconnection is present (cf. Sect. [T]i. 

However, because of the gauge invariance requirement 
which imposes that must be a flux surface, magnetic helicity 
appears as a quantity of limited practical use. In most studied 
systems, magnetic held is threading the surface d'V and the con- 
dition (B ■ dS) = 0 is not fulfilled. In their seminal paper, 
iBerger & Field! (Il984l) have introduced the concept of relative 
magnetic helicity: a gauge-invariant quantity which preserve es¬ 
sential properties of magnetic helicity while allowing a non null 
normal component of the held B through the surface of the stud¬ 
ied domain. 


2.2. Relative Magnetic helicity 

In their initial work. lBerger & Fieldid 19841) gave a first definition 
of the relative magnetic helicity as the difference between the 
helicity of the studied held B and the helicity of a reference held 
B() having the same distribution than B of the normal component 
along the surface: ((Bq - B) ■ dS) Is^y - 0. 

While the definition allows for any held to be used as the 
reference held, the potential field Bp is frequently used as a ref¬ 
erence field in the literature. Since V x Bp = 0, the potential field 
can be derived from a scalar function (p: 

Bp = V0 (8) 

where the scalar potential (p is the solution of the Laplace equa¬ 
tion A(p - 0 derived from V ■ Bp = 0. Given the distribution of 


(Bp-dS)l^n' =(B-dS)bn/ (9) 

Under the above assumptions, the potential field has the low¬ 
est possi ble energy for the given distribution of B on d'V (e.g. 
Eq.(2) of IValori et aHl201^ . In the following we will also use 
the potential field as our reference field. 

A second gauge-indepe ndent definition fo r relativ e magnetic 
helicity has been given by iFinn & AntonsenI (Il985h . definition 
used from here on in this article: 

{(A + Ap)-(B-Bp)dy. (10) 

J'V 

with Ap the potential vector of the potential field Bp = V x Ap. 
Not only is H gauge invariant, but the gauges of A and Ap are in¬ 
dependent of each others, i.e., for any set of sufficiently-regular 
scalar functions {ip\ ipp), H will be unchanged by the gauges 
transformation (A; Ap) —» (A -i- Vt/r; Ap -i- Vtpp). Let us note that 
Ap and (p correspond to two distinct solutions of the Helmholtz’s 
theorem, i.e. two distinct non-incompatible decompositions of 
Bp. 

The relative helicity in Eq. (fTOl i can first be decomposed in a 
contribution due only to B, Eq. ([T]i, one only to Bp, and a mixed 


term: 



H = 

H^-Hp+ H^ix with 

(11) 

Hp = 

r Ap Bpd^ 

(12) 




^mix “ 

r(Ap B-A Bp)dT/= r (AxAp) dS 
Jd^ 

(13) 


Let us note that this decomposition is only formal in the sense 
that each term is actually gauge dependent and only their sum is 
gauge invariant. 

Relative magnetic helicity, as defined in Eq. (fTOl l. is equal to 
the difference between the helicity of the field B and the helicity 
of its potential field (H - H^(B) - H^(Bp)) only if H^ix cancels. 
Relative helicity is in general not a simple difference of helicity, 
as for relative energy. A sufficient (but not necessary) condition 
that ensures the nullity of the mixed term is that A and Ap have 
the same transverse component on the surface, i.e.: 

A X dSla-y = Ap X dS|a<y (14) 

This condition automatically enforces the condition of Eq. (Ull 
on the normal field components. However it imposes to fink the 
ch oice of the gauge of A with that of Ap. The original definition 
of iBerger & Eieldi d 19841) cor responds to a quantity wh ich is less 
general than the one given bv lEinn & AntonsM(ll985h . It is only 
gauge invariant for particular sets of transformation: (A; Ap) —> 
(A + Ap + Vipp). 

Another poss ible decompo sition of relative magnetic helicity 
from Eq. (fTOl i is (lBergeil2003h : 


H 

- H^+ 2Hpj with 

(15) 

Hi 

- r (A-Ap)-(B-Bp)dT 

(16) 


J'V 


Hpj 

= r Ap-(B-Bp)d^ 

(17) 


'V 
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where //j is the classical magnetic helicity of the non-potential, 
or current carrying, component of the magnetic field, Bj = B-Bp 
and //pj is the mutual helicity between Bp and Bj. The field Bj 
is contained within the volume 'V so it is also called the closed 
field part of B. Because of Eq. (|9]), not only H, but also both Hj 
and Hpj are independently gauge invariant. 


Including Eqs. (I20l- l22ll in Eq. (fT9]) . the variation of magnetic 
helicity can thus be decomposed as: 


dH 

"dT 


dH 

"dT 


diss 


dH 

df 


+ Evn + Efin + EaAp + 

Bp,var 


with 


(23) 


2.3. Relative magnetic helicity variation 

Assuming a fixed domain 'V, we can differentiate Eq. (Uni) in 
time in order to study the time variations of relative helicity; 

d{A + Ap) 

a. 

"V 

a<R - BA 

d^ 


dH _ r 

df 


dt 

(A + Ap). 


■ V X (A - Ap) d^ 
5(B - Bp) 


dt 

Using the Gauss divergence theorem one obtains: 

/ 5(A + Ap) 

(A-Ap)x-^ 

\ 


dH _ r 
df Jd'i 


dt 


■dS 


f(A- 

/( 


Ap) 


+ I (A -H Ap) 


g(B + Bp) 
dt 

g(B - Bp) 
dt 


d'V 


d-V 


( 18 ) 


Combining the second and third terms, one finds the following 
synthetic decomposition of the helicity variation in three terms: 


dH 

r . -^B 

— =2 

A ■ — d^ -H 

df 

dt 


r 5Bp 

-2 

Ap • ^ dT' 


U dt 


f 

Jcrv 


(A - Ap) X 


(9(A -1- Ap) 
Jt 


■dS 


(19) 


This decomposition is only formal. Indeed, as for the decompo¬ 
sition of relative helicity of Eq. (fTTT i. none of these three terms 
are independently gauge invariant and only their sum is. The 
third term can be further decomposed using the scalar potential 
(j), Eq. (|8]l, and the Gauss divergence theorem: 



Using the Earaday law and the Gauss divergence theorem we 
also obtain: 



^d^ = - r (ExA)-dS- r B EdA^ 

dt jQry Jry 


( 21 ) 


Assuming that at the boundary the evolution of the system is 
idea l: E - (-v xBlU n/ the surface flux can be written as 
(e.g. lBerger & Eieidlll984 ): 


dH 

”57 

dH 

dt 


diss 


-2 EBd^ 


Bp,var 

F\a - 
pBn - 
Faap — 




2 I -^V-Apd'V 


-2 f (B A)v dS 
JiTV 

2 f (vA)B dS 
JdV 


5(A + Ap) 




■dS 


-2 f ^ 

JdV dt 


Ap -dS 


(24) 

(25) 

(26) 

(27) 

(28) 
(29) 


The d///df|diss term is a volume term which corresponds to the 
actual dissipation of magnetic helicity of the studied magnetic 
field (Eq.|5ll. The d///df|Bp,var term, which, despite being a vol¬ 
ume term, is not a dissipation, traces a change in the helicity 
of the potential field. As B is evolving, its distribution at the 
boundary implies a changing Bp (Eq.|3- The helicity of the po¬ 
tential field, not necessarily null, is therefore evolving in time. 
More precisely, the potential field is defined only in terms of its 
boundary values. However, this is not true in general for the 
vector potential of the potential field, because of the gauge free¬ 
dom. Hence, in general, the helicity of the potential field cannot 
be expressed as a function of boundary values only, except for 
the particular case of a vector potential without sources or sinks 
in ’V, i.e., when the Coulomb gauge is used. Therefore, the time 
variation of the helicity of the potential field necessarily contains 
both volume and flux contributions. 

All the other terms, are flux terms that correspond to the 
transfer of helicity through the surface boundary d'V. The 
Evn and Een are sometimes called the "emergence" and "shear" 
terms, but such a characterization can be misleading as their con¬ 
tributions depend on the gauge selected for A. The Eaap term is 
related to a cross contribution of A and Ap. Einally corre¬ 
spond to a flux of the helicity of the potential held. 

The d///df|diss term is the only term of the decomposition 
which is gauge invariant. All the other terms are not indepen¬ 
dently gauge invariant. It means that the relative intensity of 
these terms will be different for different gauges. Combined, 
they produce the same, gauge-invariant value of dHjdt. We will 
study the dependence of the above decomposition on the chosen 
gauge in Sect. |5] We should also note that the total flux Etot of 
relative helicity. 


Etot - Evn + Esn + EaAp + E^ , 


(30) 


- r (E X A) ■ dS = - r (B ■ A)v ■ dS -H r (v ■ A)B ■ dS (22) 
Jav JdV Jd'V 

Note that if the evolution of the system was not ideal at the 
boundary, an additional flux term depending on Enon ideal x A, 
could be added, with Enon ideal the non ideal part of the electric 
held. In the study presented here this term will de facto be esti¬ 
mated but assumed to be measured as a volume dissipation term. 


is only gauge invariant for A but not for Ap. 

Moreover, unlike d///df|diss, d///df|Bp,var is not a priori null in 
ideal MHD. It implies that dHjdt cannot be written in a classical 
conservative form since d///df|Bp,var cannot be strictly written as 
a flux term. Therefore relative magnetic helicity is not a priori a 
conserved quantity of MHD in the classical sense: its variation 
in 'V may not solely come from a flux of relative helicity through 
the boundary. Relative magnetic helicity may not be conserved 
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even if d///df|diss is small. The conservation of relative helic¬ 
ity and the dissipation of magnetic helicity are thus two distinct 
problems: the relative intensity of these terms will vary, depend¬ 
ing on which gauge is employed as we will illustrate in Sect. |5] 

2 . 4 . Relative Magnetic helicity variation with specific gauge 
conditions 

While the variation of magnetic helicity can be generally de¬ 
scribed by Eq. (l2?t for any gauge, the choice of some specific 
additional constraint on the gauge allows to simplify the expres¬ 
sion of d///df and possibly its computation. 


2.4.1. The A|fl.y = Ap|a^ condition 

We note that with the specific condition 

A|a^ = Aplfl^ , (31) 

the condition of Eq. (fT4l) is necessarily satisfied thus Eaap = 0 
and the terms Evn and Fsn can be expressed only in terms of Ap. 
Thus, the helicity variation, Eq. (|2^ . simplifies as: 


dH 

IF 



-2 f 1 
-2 f ( 

- 2 r # 

jifv dt 


E ■ B d^ -H 2 


r ^ 

dt 

£ 


V-Apd'V 


(B ■ Ap)v ■ dS H-2 (vAp)B-dS 


Ap ■ dS 


(32) 


We note that in this derivation the vector potential A is ab¬ 
sent. The condition (OTl l allows to get rid of the need of com¬ 
puting A to estimate the helicity variations. However, to derive 
H from Eq. (fTOl i both A and Ap must be computed with gauges 
coupled with Eq. (EB- Then, one must strictly control that this 
condition is enforced all over the surface of the studied system. 
This can actually be numerically challenging. 

In the present manuscript, in order to determine the helicity 
dissipation, we will compare time integrated helicity flux with 
direct helicity measurements. As our numerical method does 
not allow us to enforce Eq. OTI) . we cannot use the simplified 
Eq. (l32l l to compute the helicity variation. 


2.4.2. Coulomb gauge for potential field; V ■ Ap = 0 

If to determine the potential field one decides to choose the 
Coulomb gauge. 


V ■ Ap = 0 , 


(33) 


then there is no volume variation of the helicity of the potential 
field d///df|Bp,var = 0. The helicity variation can be reduced to 
the simplified form; 


dH 

dt 


Cond. i33] 


= - 2/1 


B-Ed^ + Ft, 


(34) 


Using the coulomb gauge for the potential field one observes 
that the variation of the relative magnetic helicity are given by a 
flux of helicity through the boundary and the dissipation term 
dH/ dtjdiss. Relative magnetic helicity to a reference field ex¬ 
pressed in the Coulomb gauge can therefore be written with a 
classical conservative equation. 


2.4.3. Boundary null Coulomb gauge 

The condition (l3^ does not enforce a unique solution for Ap. It 
is possible to further constrain the Coulomb gauge if the vector 
potential Ap satisfies the additional boundary condition; 

Ap ■ dSl^n' = 0 . (35) 

With this further constraint the flux of helicity of the poten¬ 
tial field is null, F^ = 0. Then, with conditions (l3jl l and (l3^ . the 
helicity variation thus reduces to the form; 


dH 

dt 


Cond.i. 


-if E-Bd^+ f (2(vA)B 
J'V Jd'V \ 

5(A-i- Ap)\ 

- 2(B ■ A)v + (A- Ap) X ———- ■ dS (36) 

dt i 


2.4.4. Simplified helicity flux 

Including the three conditions of previous sub-sections, i.e .,: 


i A\ay - Aplfl^ 

V-Ap =0 
Ap-dSlan^ = 0. 


(37) 


we obtai n the well known expression for the sim plified helicity 
flux (e.g. [Berger & Eieldlll984l : lPariat et al.ll2005h : 


dH 

"dT 


Cond. 1371 


-2 r E-Bd'V 

J'V 

+ 2 f ((V ■ Ap)B - (B ■ Ap)v) ■ dS 
Jd'V 


(38) 


With this set of conditions, the flux terms Fyn and Fsn are 
fully fixed. However, one should remind that these terms remain 
gauge dependent; using a gauge where the conditions of Eq. dJTl l 
are not fully enforced would lead to a different Ap and A, conse¬ 
quently, to a different distribution of helicity flux between Fyn, 
Fsn and other terms. It is therefore theoretically incorrect to 
study them independently. 

Equation (l38T l is the cl assical formulation for the helicity flux 
that has been derived by iBerger & Eieldl (Il984ll for an infinite 
plane. However in the case of a 3D cubic domain this formula¬ 
tion is only valid if all t he conditions of Eq. (IJTI) are satisfied. 
While l^ng et al.l (1201 3h have used Eq. dMt to compute the he¬ 
licity flux, it remains to be determined if all conditions (iJTl i are 
fulfilled when computing the volume helicity. Their relatively 
high level of non-conservation (3%) in the ideal phase of the 
evolution of their system may be related to this discrepancy. Fi¬ 
nally, while conditions of Eq. dJTl i can drastically simplify the 
estimation of the helicity variation, they strongly constraint the 
numerical implementation of A and Ap. Fast, precise and prac¬ 
tical numerical computation of the vector potentials may require 
a different choice a gauge. 


3. Methodology 

The present section describes the methodology employed in our 
numerical experiments. In Sect. 13.11 we present the general 
method used to estimate the relative magnetic helicity conser¬ 
vation and the magnetic helicity dissipation. Then in Sect. 13.21 
we describe how we practically compute the volume helicity and 
its flux. Finally, in Sect. 13.31 we present the numerical data set 
considered for the helicity conservation tests. 
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Fig. 1. Snapshots of the magnetic field evolution during the generation of the jet. The red field lines are initially closed. The green and white 
field lines are initially open. All the field lines are plotted from fixed footpoints. The red and white field lines are regularly plotted along a circle of 
constant radius while the green field lines are plotted along the x axis. At ? = 900 the system is in its pre-eruption stage. It is close to the maximum 
of energy and helicity. All the helicity is stored in the close domain. At f = 1000 the system is erupting. Numerous field lines have changed of 
connectivity as one observe open red field lines and closed white field lines. Helicity is ejected upward along newly open reconnected field lines. 
At f = 1500 the system is slowly relaxing to its final stage. 


3.1. Estimators of magnetic heiicity conservation 

Following lYang et al.l (l2013h . we will compute the volume vari¬ 
ation along with the flux of magnetic helicity. From two suc¬ 
cessive outputs of the studied MFID system, corresponding to 
two instant t and t', separated by a time interval At we will di¬ 
rectly compute their respective total helicity and 

in the volume 'V using the method of Sect. I3.2I and then, the 
helicity variation rate between these 2 instants: AH^jAt - 
{Hry(T') - H<v{t))IA t. We will simultaneously estimate the dif¬ 
ferent sources of fluxes F# through the surface of the domain, 
with # the different contribution to the total flux in Eq. (l30l) . 

We will also time integrate the helicity fluxes at the boundary 
Hav- 


Hgv(t) - 

f Ftotir) dr 

Jo 

(39) 

Hd'V.dt) - 

r F#(T)dT 

Jo 

(40) 


with # the different contribution to the total flux of Eq. (l30l) . 

One should note that two very different methods are used 
to compute H<v and Hgry. To derive (and AHryjAt) from 
Eq. (ITOt . only three components of the magnetic held B in the 
whole domain 'V are needed: A,Bp, Ap are derived from B (cf 
Sect. I3.2k Eor Hgry (and Ftot), only data along the boundary d'V 
are required. Eurthermore, the helicity flux estimation requires 
the knowledge of the three components of the velocity held on 
d'V (to compute Fvn and Fbi,)- These quantities are not used to 
compute H<v. The methods of estimations of H<v and Flgy are 
thus fully independent. 

A physical quantity is classically said to be conserved when 
its time variation in a given domain is equal to its flux through 
the boundary of the domain (cf. Sect. [T]l- To study the twofold 
problem of the conservation of relative magnetic helicity and the 
dissipation of magnetic helicity, we will use two quantities, Cr 
and Cm, respectively. 

Eor relative magnetic helicity to be perfectly conserved one 
should have Hry - Hgry^ or dHIdt must be equal to Ftot, i.e. the 
volume terms of the relative helicity variation should be null. By 
estimating, Cr, 


Cr = 


AH^ 

At 


dH 

dt diss 


dH 


df 


Bp,var 


(41) 


we can determine the level of conservation of relative magnetic 
helicity. As already noted in Sect. 12.31 even if B evolves within 


the ideal MHD paradigm, the relative magnetic helicity evolu¬ 
tion cannot be written with a classical equation of conservation 
as the term in d///df|Bp,var is not a flux integral through a bound¬ 
ary and generally does not vanish. 

Secondly, we aim to determine the dissipation of the mag¬ 
netic helicity of the studied held (Ea.l24l). i.e. by estimating Cm 
equal to: 


r - 


AHry 

At 


~ Ffnt ~ 


AHry 


At 


Bp,var 


dH 

dt 


diss 


(42) 


Our estimation of Cm is done independently of the estimation of 
dHldt\aiss- Our method thus does not require the knowledge of 
the electric held E, which is a secondary quantity in most MFID 
problems. The dissipation is thus estimated in a way which is 
completely independent from the non-ideal process, e.g. the pre¬ 
cise way magnetic reconnection is developing. 

As the estimators measure both physical helicity variations 
and numerical errors, we will use different non-dimensional cri- 
terions to quantify the level of helicity conservation and the pre¬ 
cision of our measurements. We will compute the relative accu¬ 
mulated helicity difference en: 


enCO - 


Hyt) - Hsrv(t) 
Href 


X=o 

Href 


(43) 


with Href a normalizing reference helicity, physically significant 
for the studied system (e.g. the maximum H^y value in the stud¬ 
ied interval). 

At each instant. Cm expresses the rate of dissipation of he¬ 
licity, i.e., our numerical estimation of Eq. (l24ll . It can include 
artefact fluctuations due to the numerical precision of the flux 
estimation and the time derivation of H<y. In addition, as helic¬ 
ity is a signed quantity both positive and negative helicity can 
be generated by non-ideal effects. It may be relevant to know 
the time integrated absolute variation of helicity, as a function of 
time. Hence, we define another metric, ecm, 


^Cm(0 ~ 


!U\Cn.\dr 

Jl^lAHey/Atldr 


(44) 


where the absolute values guarantee that we take an upper limit 
of the dissipation. Along with en, fcm allows us to evaluate the 
level of dissipation and the numerical precision of our measure¬ 
ments. The measure of low value of ecm and en can thus pro¬ 
vide a clear demonstration of the level of conservation of relative 
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magnetic helicity. For two methods having similar eu, a higher 
value of ecm would indicate measurements presenting larger nu¬ 
merical errors. 




Fig. 2. Top panel: Time evolution of the different magnetic ener¬ 
gies relative to their respective initial values: total (SE^y, black dashed 
line), of the potential magnetic fields (SE-y^p, purple line), of the non 
solenoidal component (6E^i,„ orange line), and free magnetic energy 
(SEryj, red line). Bottom panel: Relative magnetic helicity {Hry, black 
dashed line) and its decompositions, Eqs. ( fTTTTsI ). computed with the 
practical DeVore method: proper helicity blue line), potential field 
helicity (Hp, purple line), and the mixed helicity (//mix> green line) ; 
current-carrying helicity (H^, red line), mutual helicity between poten¬ 
tial and current-carrying fields (//pj, orange line). In both panels, as in 
all figures hereafter, the dot-dashed vertical line at f = 920 indicates the 
transition between the quasi-ideal-MHD/energy accumulation and the 
non-ideal/jet generation phases. 


3.2. Volume helicity and flux computations in the DeVore 
gauge 

In order to compute magnetic hel icity in the v olume , Hry, we 
will use the method presented in IValori et ID ( I 2 OI 2 I) . All the 
vector potentia ls will be computed using the gauge presented in 
iDeVorel (I 2 OOOI) . with no vertical component: 

=0 (45) 

Under this assumption, the vector potentials can be com¬ 
puted i n the volume usin g a ID integral (c.f. Equations (10) or 

(11) of IValori et alJl201^ with a 2D partial differential equation 
to be s olve at the bottom or top boundary (c.f. Equations (9) or 

(12) of IValori et al.ll20f^ . We tested that we did not obtained 
significant differences on the helicity conservation and dissipa¬ 
tion properties whether the integration was done from the top or 
the bottom boundary of the system. In the present manuscript 
the results were obtained with the integration performed from 
the top boundary. 

While other gauge choices could be explored, choosing the 
DeVore gauge, Eq. (l45l l. is here motivated by the fact that it is nu¬ 
merically efficient and convenient. This selection leaves a free¬ 
dom for the gauge of A and Ap which can be independent, i.e. 
not linked as in Eq. OTT i. In our computation of Ap and A, we 
use gauge freedom to additionall y enforce that the two vectors 
are equal at the top (see Eq. 29 of lValori et al.ll201^ : 

A(ztop) ~ Ap(ztop) (46) 

We still have a freedom on An(ztnn) as expressed by the 2D partial 
differential equation (20) of lValori et al.l (l2012l) . Here, we select 
their particular solution expressed by their equations (24,25). 
With this additional choice, both Ap and A are uniquely defined 
(modulo a constant) by vertical integration starting from the top 
boundary. 

Practically, we first determine the potential fiel d by solving 
the so lution of the Laplace equation (17) for (p of lValori et al.l 
(l2012h . Then we can compute the potential vectors from a direct 
ID integr ation of starting from the top boundary (c.f. Sec¬ 
tion 3.3 of lValori et al.ll20O . We refer to this method below as 
the "practical DeVore method" since it is efficient and easy to 
implement. 

While the condition (l46T l holds at the top boun dary it does 
not ho ld at the other boundaries. Indeed, Eq. 13 in IValori et al.l 
(I 2 OI 2 I) can be used to show that the difference between the yalue 
of Ap and A at the bottom boundary is equal to z x '(B - 
Bp)d'T'. It is important to note that the conditions of //mix = 0 
and Ap ■ dSlaq/ = 0 are neyer enforced in this gauge. The relatiye 
helicity terms //mix and //p and the helicity flux term /^aap and 
Fit, can thus neye r be considered null (see also Section 3.4 in 
IValori et^l2012i for more details). Therefore, for the practical 
DeVore method, the helicity yariation of the system is giyen by 
the general formula of Eq. (|2^ . 

Alternatiyely, we will determine the helicity assuming a 
Coulomb gauge for Ap only. The DeVore and Coulomb gauge 
are indeed compatible for a potent ial field. We can so lye the 2D 
partial differential equation (20) of lValori et al.l (12012h as a Pois¬ 
son problem to obtain Ap (see their Eq. (41), but translated to 
the top boundary). This implies by construction that Ap is si¬ 
multaneously respecting the DeVore, Eq. (1451) . and the Coulomb 
gauge, Eq. (l3^ . In the following, we will refer to this method as 
the “DeVore-Coulomb method”. 

While still following the DeVore condition (|45]) . A computed 
in the DeVore-Coulomb method will be different from A com¬ 
puted in the practical DeVore method because the boundary con¬ 
dition, Eq. (l46l l is different. In particular, the distributions of A 
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at the bottom boundary will be significantly different with each 
methods, hence leading to very different values of F^n and Fvn- 
In AppendixiAl we present an additional test, with Ap computed 
in the Coulomb gauge but where A is not satisfying the condition 

(E 6 ll. 


3.3. Test data set 


In order to test the helicity conservation, we employ a test-case 
3D M HD numerica l simul ation of the generation of a solar coro¬ 
nal iet IPariat et alJ (l2009l) . Figure [T] presents snapshots of the 
evolution of the magnetic field. The simulation assumes an ini¬ 
tial uniform coronal plasma with an axisymmetric null point 
magnetic configuration (left panel). The magnetic null point 
is created by embedding a vertical dipole below the simulation 
domain, and adding an uniform volume vertical magnetic field 
of opposite direction in the domain. The null point present a 
fan/spine topology, dividing the volume in two domains of con¬ 
nectivity, one closed surrounding the central magnetic polarity, 
and one open. 

The computation are performed with non-dimensional units. 
We analyze the time evolution of the magnetic field from f = 0 
to f = 1600. The time steps between two outputs is At = 50 for 
t e [0;700], during the accumulation phase, and. At = 10 for 
t 6 [700; 1600] during the dynamic phase of the jet. The anal¬ 
ysis volume 'V is a subdomain of the larger discretized volume 
employed in the original MHD simulation: [-6,6] in x-, [-6,6] 
in y-, and [0,12] in the z-directions, thus only tak ing into ac¬ 
count the region of higher resolution (c.f. Figure 1 of lPariat et al.l 
120091) . 


The ideal MHD equations are solved with the AR MS code 
based on Flux Corrected Transport algorithms (FCT, iDeVorel 
1199 Ih . The paralle lisation of the code is ensured with the 
PARAMESH toolkit ( MacNeice et al.l2000h . In the original sim¬ 
ulation of lPariat et al. ( 20091) . reconnection is strongly localized 
thanks to the use of adaptive mesh refinement methods at the lo¬ 
cation of the formation of th e thin current sheets involved in re¬ 
connection (c.f. Appendix of lKamen et al.ll2()T^ . In the present 
paper in order to keep the resolution of the domain constant we 
have however switched off adapti vity, so the resolu tion is equal 
to the initial one (as in Figure 1 of lPariat et al.l200^ . throughout 
the whole simulation. 

The variation of energies relatively to their initial value, 
6Erv{t) - Ery{t)-Ery{t - 0) is displayed in Fig.|2l top panel. The 
index 'V indicates that the energy is computed by a volume inte¬ 
gration. The total energy, E^y, in the domain can be decomposed 
as: 


E - + E^^ j + , (47) 

where E<v.p and E<v, j are the energies associated to the poten¬ 
tial and current-carrying solenoidal contributions and Ery „s is 
the sum of the n onsolenoidal contributions (see Eqs. (7,8) in 
IValori et'^l2013L for the corresponding expressions). Initially 
the system is fully potential and E<v - E^y p. Energy is injected 
in the system by line-tied twisting motions of the central po¬ 
larity. The axisymmetric boundary motions are preserving the 
distribution of at the bottom boundary. Magnetic free energy 
and helicity accumulates monotonically increasing the twist in 
the closed domain (Eig.[Tl central left panel). The potential field 
energy E<y^p decreases slightly during the accumulation phase 
because of the bulge of the central domain which is changing 
the distribution of the field on the side and top boundaries. The 
E'y^ns term is almost constantly null, an indication of the excel¬ 
lent solenoidality of the system 


Eventually, around t ^ 920, the system becomes vio¬ 
lently unstable: magnetic reconnection sets in and the closed 
twisted field lines reconnect with the outer open field lines. 
A steep decrease of the free magnetic energy is observed in 
Fig.m 83% of the maximum free magnetic energy is dissi¬ 
pated/ejected/transformed. Through reconnection, twist and he¬ 
licity are expelled from the central domain, inducing a large 
scale kink wave which exits through the top boundary (Eig. [1] 
central right panel). This large scale non-linear magnetic wave, 
simultaneously compresses and heats the plasma, inducing the 
generation of an untwistin g jet that can observationally be in¬ 
terpre ted as a blowout jet (iPatsourakos et aljl2008t iPariat et ^ 
l2015l) . The driving motions had been slowly ramped down at 
the time of the trigger of the jet so that few energy and helicity 
are inje cted from the lowe r boundary after the jet onset (cf. Eig- 
ure 6 of lPariat et al.ll20d^ . In the final stage the system slowly 
relaxes to a configuration similar to its initial state, with the po¬ 
tential field energy being similar to its initial value, and only few 
field lines remaining twisted, next to the inversion line (Pig. [T] 
right panel). 

This simulation thus presents two distinct phases typical of 
active events: before t ^ 920, a phase with a slow ideal ac¬ 
cumulation of magnetic helicity and energy, and after t ^ 920, 
an eruptive phase of fast energy dissipation and helicity trans¬ 
fer involving non-ideal effects. In the first phase, the system 
behaves very close to ideality a s demonstrated in a be nchmark 
with a strictly i deal simulation ([Rachmeler et alJl2010h . In the 
non-ideal phase. IPariat et alJ (l2009l) showed that 90% of the he¬ 
licity w as eventually eject ed through the top boundary by the 
jet, and IPariat et M (1201 Ol) showed the high reconnection rate 
processing the magnetic flux during the jet. These two phases 
allow us to test the conservation of helicity in two very distinct 
paradigms of MHD. 

In the following we test the helicity conservation with the 
above MHD simulation. We first use the practical DeVore 
method (see Sect. Si. We next test the effect of the gauge 
choice on our results, by using the DeVore-Coulomb method 
(see Sect. |5]l. We will show how the gauge choice possibly af¬ 
fects the evolution of each terms. 


4. Magnetic helicity conservation in the practical 
DeVore gauge 

Our first study of the helicity variations is performed in the prac¬ 
tical DeVore gauge. The only assumptions on the potential field 
and the vector potentials are given by Equations (0, (1451) & (l46l l. 
This corresponds to a general case where the helicity variation 
of the system is provided by Eq. (|2^ . 


4.1. Helicity evolution 

Figure |2] bottom panel, presents the evolution of the relative 
magnetic helicity H<y in the system. Similarly to magnetic en¬ 
ergy the two phases of the evolution are clearly marked. The first 
phase cotTesponds to a steady accumulation of magnetic helic¬ 
ity, while the second corresponds to the blowout jet. The latter 
is associated with a steep decrease of magnetic helicity. As the 
system relaxes the magnetic helicity value oscillates. These os¬ 
cillations are related to the presence of a large scale Alfvenic 
wave which is slowly damped after the jet. This oscillations can 
also be seen in the total energy but with a much smaller relative 
amplitude (so small that E^ still decreases monotonically). 
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Fig. 4. Total helicity flux, Ftot lEg. I30t . and the the terms compositing it, Eqs. I l26l -I29t, through the different boundaries computed in the 
practical DeVore Gauge, (a): Ftot; (b): (c): Faap; (d): Fenl (e): Fvn. In each plot the dark line correspond to the sum of the flux through all 

the boundaries while the color lines correspond to a flux through a particular boundary (purple and blue: left and right x sides, green and orange: 
front and backy sides; yellow and red: bottom and top z sides). 


The right panel of Fig. |2] also presents the decomposition 
of the relative magnetic helicity, H<v, in f/mix and of 
Eqs. (ED-inil . During the accumulation phase, H<v is domi¬ 
nated by i/tnix while is almost constantly null. During the 
jet, one can see strong fluctuations of the relative importance of 
these terms. //„, //niix and eventually present contributions of 
similar amplitude. However, because these terms are not gauge 
invariant, their value in a different gauge might be quite different 
(cf. Sect.|5]). 

The decomposition of H<v with //j and //pj of Eqs. (na-inii 
is also plotted in Eig. |2l bottom panel, //j has an evolution com¬ 
parable to 6E<v y It captures most of the helicity evolution both 
during the accumulation and jet phases. In contrast, the mu¬ 
tual helicity between the potential and the current carrying fields 
contains mostly oscillations. Therefore, Hy which is a gauge 
invariant quantity, is a promising quantity to analyze during a 
jet/eruption. 

4.2. Helicity fluxes 

We compute the time variations AHr^jAt of the helicity deter¬ 
mined with the volume integration method, and compare it with 
the different terms of the relative magnetic helicity flux through 
the whole system boundary (Eig.[3]|. This shows that the helic¬ 
ity variation AH<vlAt is very closely matching the curve of Ftot, 
indicating that the variation of helicity in the domain are tightly 


related to the flux of relative helicity through the boundary. The 
core results of this study is that indeed magnetic helicity is very 
well conserved in the studied simulation, both during the quasi- 
ideal-MHD and non-ideal phases. 

During the non-ideal/jet phase, strong fluctuations are ob¬ 
served for all terms but the Eaap term. The later is constantly 
negligible relatively to the others. On the other hand we see that 
the term displays important fluctuations, frequently of similar 
amplitude and opposite sign to E’en- The term clearly cannot 
be neglected in this gauge. 

The analysis of the flux of each terms through each individ¬ 
ual boundary is insightful (Eig.|4]i, though one must bear in mind 
that the plotted terms are not gauge invariants. Although its am¬ 
plitude is extremely small, a finite flux of Eaap is present only 
at the bottom boundary during the whole evolution of the sys¬ 
tem; there is no flux on the sides because of the DeVore gauge 
(Ea.l45]). and no flux on the top because of the imposed condition 
of Eq. (Hhl) . 

During the ideal phase, the flux of helicity is completely 
dominated by the Esn term originating from the bottom photo- 
spheric boundary. This is consistent with the fact that the system 
is indeed driven by horizontal shearing motions at the bottom 
boundary. No remarkable helicity flux is observed at the other 
boundaries during this period. Helicity is thus accumulating in 
the volume W'. 
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Fig. 5. Difference between the helicity 
variation rate and the boundary helicity 
flux, computed in the practical DeVore 
gauge. The relative helicity conserva¬ 
tion criterion Q (red line, Eq. ED and 
the magnetic helicity dissipation crite¬ 
rion Cn, (blue line, Eq. I42t are plotted 
relatively to the helicity variation rate 
(AH^/At, black dashed line, left panel) 
and in their own amplitude range (right 
panel). The insert in the right panel 
present the difference between Cr and 
Cm, i.e. the potential helicity volume 
variations d///dfp,var- 



Fig. 3. Comparison of the helicity variation rate and the helicity flux 
computed in the practical DeVore gauge. The helicity variation rate 
(AH^vlAt, dashed line) is derived from the volume integration method. 
The total helicity fluxes through the whole surface, (TV, of the domain 
is Ftot (red line). This flux is composed of, Eqs. i26l - l29t : Fvn (green 
line), Fbh (blue line), (orange line), Faap (purple line). 

During and after the jet (f > 920) important helicity fluxes 
are noted in the side and top boundaries, while the bottom flux 
is now negligible as the boundary flows have been ramped down 
in amplitude. A large flux of helicity occurs at the top boundary 
(red curves), dominated by the Fan term and to a lower extend by 
the Fvn term. This corresponds to the ejection of he licity by the 
jet, driyen by a large-scale non-linear torsional waye (iPariat et al.l 
(20091 l2015h . Fvn is peaked at the time of the passage of the 
bulk of the jet through the top boundary. The flux of Fb,, and 
Fvn through the side boundaries, while present, is comparatiyely 
small. Howeyer, the side boundaries see the transit of important 
flux F^. No specific side boundary is dominating the total yalue 
of F^. Due to the DeVore gauge (Eq. l45T l. F^ is null at the bottom 
and top boundaries. 


We conclude that, computed with the particular DeVore 
method, the total flux of helicity Ftot during the jet consists of 
complex transfer of helicity through all the boundaries of the 
system. 

4.3. Helicity conservation 

In order to better estimate the helicity conseryation and dissipa¬ 
tion, we haye plotted in Fig.|5]the criterions Cr and Cm (see their 
definition in Sect. 13.II) . We obserye that these two criterions are 
almost equal. They differ only by the term describing the yol- 
ume yariation of the potential helicity d///df|Bp,var of Eq. (l26l l. 
Our calculation finds that d///df|Bp,var is extremely small com¬ 
pared to the yariation of Cr and Cm (see inset in the right panel 
ofEig.H), eyen though we are not explicitly enforcing V -Ap = 0. 
This term is smaller than the combined effect of the real helic¬ 
ity dissipation d///df|diss and the numerical errors on the yolume 
and flux helicity measurements. 

The curye of Cr demonstrates that magnetic helicity is ex¬ 
tremely well conseryed during the ideal-MHD phase of the sim¬ 
ulation and that it is also yery well conseryed during the non¬ 
ideal phase. Moreoyer, during the ideal phase the amplitude of 
|Cml does not exceed 0.029 which is 1% of the maximum ampli¬ 
tude of helicity yariation during the period. At the end of this 
period, one also haye ecm(f = 920) <1%, t hus helicity is yery 
weakly dissipated, as theoretically expected (IWoltieilll958h . 

During the non-ideal phase, Cr and Cm show high frequency 
fluctuations around the null yalue, decorrelated from the fluc¬ 
tuation of helicity in the system. Our analysis indicates that 
while these fluctuations could originate from the real physical 
term d///df|diss, they are in fact dominated by the numerical pre¬ 
cision on the estimation of Ftot- From the same simulations data, 
we will see in Sect. |5] that the computation with the DeVore- 
Coulomb method reduces these fluctuations. 

In Fig. |6] we haye plotted the yariation of magnetic energy 
and helicity in the system computed with a yolume integra¬ 
tion and from the integration of the Poynting and helicity fluxes 
through the boundaries of the system. During the ideal-MHD 
phase both magnetic energy and helicity are well conseryed, 
their yolume yariations being equal to their boundary fluxes. 
During the non-ideal phase, while magnetic helicity is still yery 
well conseryed magnetic energy is clearly not. When the jet is 
generated, the magnetic energy quickly decreases: part of it is 
ejected through the top boundary by the jet, but for most part 
is dissipated in the reconnection current sheet and transformed 
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Fig. 6. Left panel: Comparison of 
the evolution of magnetic energy (6E<y, 
black line) in the volume with the time 
integration of the Poynting flux through 
the whole surface of the domain (6Egr^;, 
red line). Right panel: Comparison 
of the evolution of helicity obtained by 
volume integration {H<y, black dashed 
line) with the time integration of the 
helicity flux through the whole surface 
of the domain (Hs'v, red line), com¬ 
puted with the practical DeVore gauge. 
Their difference is plotted in orange on 
a different range of amplitude (c.f. right 
axis). 


in other forms of energy. When the simulation is stopped, we 
determined that about 17% of the magnetic energy injected in 
the system by the bottom boundary motions remains in the sys¬ 
tem, 21% is directly ejected with the jet and the 62% remain- 
ing are dissipated /transformed in other form of energy (see also 
iPariat et ^ |2009|) As expected during reconnection, magnetic 
energy is strongly non-conserved. 

The situation for magnetic helicity is very different. At 
the end of the ideal phase the maximum difference between 
//-y and is of 3.5 units. During this period a total of 
Href - max(H) — 1265 units have been injected and only a 
fraction e^it = 920) = 0.3% is lost. The ideality of the system 
is thus very well maintained. During the non-ideal phase, the 
maximum difference between //y and Hg^y is equal to 58 units 
(orange line, Fig.|6] right panel) and enit — 1600) = 4.5%. The 
jet is able to transport away a huge fraction, 77% of the helicity 
of the system. Finally, about 19% of the helicity remains in the 
system at f = 1600 while H<v is still oscillating slightly (Fig.|6] 
right panel). 

Summarizing, with a generi c gauge , for a solar-like active 
events, we have thus confirmed iTavloj (Il974t) hypothesis that 
magnetic helicity is very well conserved even when non-ideal 
processes are acting. The relative helicity dissipation is 15 times 
smaller that the relative magnetic energy dissipation. How¬ 
ever, with the practical DeVore Gauge, we observed that Cm 
is strongly fluctuating, possibly limited by numerical precision. 
We now want to test whether we can obtain better results with 
a different gauge, simultaneously allowing us to test the gauge- 
invariance of the different helicity variation terms. 


5. Magnetic helicity conservation in the 
DeVore-Coulomb case 

Relative magnetic helicity has been defined as a gauge-invariant 
quantity (see Sect. 12.31) . However, the surface flux of relative 
helicity is not, neither are the individual terms that defines the 
flux. In the following we will study the influence of computing 
magnetic helicity using a different gauge. The vector potential of 
the potential field, Ap, is now computed in the DeVore-Coulomb 
gauge, but not A (since the Coulomb gauge is only compatible 
with DeVore gauge for a potential held). Because we also im¬ 
pose the condition of Eq. (1461) . A is however also recomputed 
during our DeVore-Coulomb method (see Sect. 13.21 ). The case 
where Eq. (l46l ) is not enforced is briehy discussed in Sect. |A] 


5. 1. Gauge dependance 

In Eig. |7] we have plotted the evolution of the relative magnetic 
helicity, H<v in the DeVore-Coulomb gauge for Ap. The left 
panel shows the comparison with the same quantity computed 
using the practical DeVore gauge, used in the previous section. 
We observe that the two curves matches almost perfectly. Their 
maximum differences is at most 3.5 units, and their maximum 
relative helicity difference is of 0.3%. The gauge invariance is 
thus very well respected for our estimation of the relative mag¬ 
netic helicity. 

In Eig.|7] right panel, the decomposition of the relative mag¬ 
netic helicity is plotted. Unlike with the gauges of the practical 
DeVore method (cf. Eig. |2]), is now almost constantly null. 
It is not strictly null since we are not imposing Ap.dSla^ = 0 on 
the side boundaries. We see that when the jet is developing 
fluctuates slightly (barely visible on Eig.|7]). The helicity is thus 
distributed differently between //mix and in the DeVore- 
Coulomb case compared to the practical DeVore gauge. As ex¬ 
pected these helicity terms are indeed not gauge invariant. Most 
of the helicity that was carried by for the practical DeVore 
gauge is now carried by while Hmm remains very similar in 
both gauges (compare the green curves in Eigs. |2] and |7]). In 
general, we have to expect a completely different distribution of 
//m, //mix and //p in other gauges. On the other hand, the quan¬ 
tity //j and //pj remains equal in both gauge computation (with 
the same precision as H<y). This is expected since //j and //pj are 
gauge-invariant quantities. 

Eor the helicity fluxes, in the DeVore-Coulomb case the time 
variations AH<ylAt of the helicity follows tightly the helicity 
flux through the boundary (Fig. [3 left panel). While AH^yjAt 
is gauge invariant, the different contribution of Ftot are not. By 
comparison with Fig. |3] one observes that /^Bn, /^Vn, /^^ have a 
very different evolution. While was presenting fluctuations 
of the same amplitude as /^tot in the practical DeVore gauge, this 
term is now very weak throughout the evolution of the system 
with the DeVore-Coulomb method. The Fen term is dominating 
Ftot in the DeVore-Coulomb case both during the ideal MHD 
phase (at the bottom boundary) and during the non-ideal period 
(at the top boundary). Fvn is weaker than with the practical De¬ 
Vore gauge. The Faap term remains negligible in both cases, 
although with a different choice of gauge, Faap can contributes 
significantly to Ftot (see Appendix lAl). 

The impact of the gauge dependance is more strikingly illus¬ 
trated when looking at the time integrated helicity fluxes //fl<y,# 
(Eq. l40l ) through the boundaries. Their evolution is presented 
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Fig. 7. Left panel: Evolution of 
the magnetic helicity computed with the 
practical DeVore method (dashed line) 
and with the DeVore-Coulomb method 
(red line). Their difference is plotted 
with an orange line {c.f. right axis). 
Right panel: Volume magnetic helicity 
(Hn/, black dashed line) and its decom¬ 
position in the DeVore-Coulomb case. 
The plotted curves are similar to Figure 
121 bottom panel. 


in Fig. |9] for the two gauges computations. While Hr^ remains 
gauge invariant (Fig.|2l left panel), as expected the i/av,# present 
different profiles in the different gauges, is small when 

computed with the DeVore-Coulomb gauge while it presents 
large amplitudes in the practical DeVore gauge. Hg^ vn and 
Hgry Bn present smaller mean absolute values when Ap follows 
the DeVore-Coulomb conditions. Their ratio is also strongly de¬ 
pendent on the set of gauges employed. 

Numerous studies have computed Hgry y^, & Hg^gn and 
followed their time evolution in observe d active regions (e.g. 
IZhang etTni2012l: iLiu et alJl2013L l2014bllah . These terms have 
been incorrectly called "emergence" and "shear" terms and phys¬ 
ical insight have been tried to be extracted from their values and 
respective ratio. However, as these terms are not gauge invari¬ 
ant, one must question the pertinence of such results. Hgryyn 
and Hgry B„ cannot be studied independently as, for a given 
V, their intensities and respective values can be simply mod- 
ified by a change of gaug e. This extends the conclusion of 
iDemoulin & Berger! (l2Q03» which showed that only the sum of 
Hgyy„ and Hgy Bn can be derived when only the tangential ve¬ 
locity components are known. More generally, the study pre¬ 
sented here shows that, even when the full velocity field on the 
boundary is known, Hgyyn and Hgy Bn can present different val¬ 
ues depending on the gauge used for the computation. In sum¬ 
mary, only the sum Hgy of all the flux terms, Hgy^#, is carrying 
a meaningful information. 

5.2. Helicity dissipation 

The estimators Cr and Cm (cf. Sect. 13.1b enable us to study the 
helicity conservation and dissipation. Since d///df|Bp,var is null 
because of the Coulomb condition on the potential field. Cm and 
Cr are equal to our numerical precision. There is no volume 
helicity variation due to the potential field. The conservation 
of the relative helicity is only limited by the actual dissipation 
helicity d///df|diss- The curve of Cr for the DeVore-Coulomb 
computation confirms that magnetic helicity is gauge-invariantly 
well conserved during the simulation (Fig.0 right panel). 

In fact the computation in the DeVore-Coulomb case even 
improves our estimation of the helicity dissipation. During 
the quasi-ideal phase, the computation of Cm in the DeVore- 
Coulomb case is equal to its value with the practical DeVore 
case. However we observe that in the non-ideal phase Cm 
presents smaller oscillations when computed in the DeVore- 
Coulomb case. Theoretically Cm, being equal to the helicity dis¬ 


sipation, d///df|(iiss, should be gauge invariant. However we ob¬ 
serve that Cm can change by a factor 2 when computed with the 
different methods. With the practical DeVore method, the fluctu¬ 
ations of Cm peak at 25% of the maximum amplitude of helicity 
variation during the non-ideal phase, while the peak is only equal 
to 4.5% of the amplitude in the DeVore-Coulomb case. The non- 
dimensional criterion, ecm at the end of the simulation is equal 
to 7% in the practical DeVore case while it is limited at 2.7% 
with the DeVore-Coulomb method (see the inset in Fig. [H right 
panel). 

While the measure of Hy is done with a high precision (see 
above), the estimation of the helicity flux, Cmt induces more nu¬ 
merical errors. The high frequency oscillations observed in the 
different terms contributing to Cmt are the symptom of the lower 
level of precision. The criterion Cm is in fact equal to the helic¬ 
ity dissipation plus the numerical errors in the volume helicity 
variations and in the helicity flux. What Cm is in fact providing 
is an upper value for the helicity dissipation. In the DeVore- 
Coulomb case. Cm is likely providing a smaller value, hence a 
better constraint on d/Z/dfldiss, because the Fxoi is dominated by 
the errors on Ffin, while in the practical DeVore case, the errors 
of the fluctuating Feh, Fy„ and will sum up with comparable 
magnitudes. 

The measure of the helicity flux through the boundary is 
the limiting factor that does not permit us to reach numerical 
precision during the non-ideal phase. It is likely that a higher 
spatial resolution sampling of the velocity field on the bound¬ 
ary sides will further improve the helicity dissipation estimation. 
Our present computation of Cm is therefore only an upper bound 
on the real helicity dissipation d///df|diss. Because of its lower 
ecm, the DeVore-Coulomb method allows us to better bound the 
helicity dissipation. 

In DeVore-Coulomb computation, during the ideal-MHD 
phase, the maximum difference between Fly and Hgy is of now 
of 4.3 units and e\i{t - 920) = 0.3%. As theoretically expected, 
the dissipation of magnetic helicity is extremely small when only 
ideal process are present. It actually demonstrates that the ide¬ 
ality of the system is very well maintained by the numerical 
scheme during that phase. At the end of the non-ideal phase, 
the maximum difference between Hy and Hgy is equal to 28.2 
units (orange line. Fig. |9] right panel) and the relative amount of 
helicity dissipated is smaller than e\i{t - 1600) = 2.2%. 

In absolute value, the dissipation of magnetic helicity is thus 
very small compare to the magnetic helicity which is ejected 
away. The dissipated helicity represents less than 3% of the he- 
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Fig. 8. Left panel: Comparison of the 
helicity variation rate and the helicity 
flux integrated through the six bound¬ 
aries of the domain in the DeVore- 
Coulomb case. The plotted curves are 
similar than in Figure Right panel: 
magnetic helicity dissipation criterion 
Cm, Eq. (|42}, plotted for the practical 
DeVore (black line) and the DeVore- 
Coulomb (red line) methods. The in¬ 
sert shows ecm, Eq- (El), with the same 
color convention. 




Fig. 9. and evolution in the 
system computed with different gauge. 
The left panel has been computed with 
the practical DeVore method while the 
right panel has been derived with the 
DeVore-Coulomb method. The helic¬ 
ity variation black dashed line) 

is derived from the volume integration 
method. The time integrated helicity 
is (red line) which can be decom¬ 
posed in the contribution of its different 
terms: (blue line), (green 

line), (orange line) and 

(purple line). The orange line in the 
inserts correspond to the difference be¬ 
tween H<v and Hgn/. 


licity carried away by the jet. If one want to track the helicity 
evolution of this jet system, the helicity dissipation only rep¬ 
resent a very minor contribution, more than one order of mag¬ 
nitude smaller than the helicity staying in the system and the 
ejected helicity. 

This low helicity dissipation is also to be compared to the rel¬ 
ative decrease of magnetic energy of 62% which is developing 
simultaneously. The reconnections generating the jets are thus 
dissipating/transforming magnetic energy > 30 times more effi¬ 
ciently than magnetic helicity is dissipated. While vast amounts 
of energy are lost, magnetic helicity is barely affected by the 
non-ideal MHD processes at play. 


6. Conclusion 

Based on the property that magnetic h elicity present s an in¬ 
verse cascade from small to large scale, iTavlod (Il974l) conjec¬ 
tured that magnetic helicity, similarly to pure ideal MHD, is 
also effectively conserved when non-ideal processes are present. 
Because of the inherent difficulties to measure magnetic helic¬ 
ity, tests of this conjecture have so far been very limited (cf. 
Sect. [p. The theoretical de velopment of relative magnetic he¬ 
licity ^BergeT&[Rei3|T983), as well as the publication of recent 
methods to actually measure relative magnetic helicity in general 
3D datasets dThalmann et al.l201 lUValori et al.l20 121 : 1 Yang et al.l 
1201 3l) are now opening direct ways to test the conservation of 
magnetic helicity. 


In this manu script, we hav e performed the first precise and 
thorough test of iTavloil (Il974l) hypothesis in a numerical simu¬ 
lation of an a ctive solar-like eve nt: the imp ulsive generation o f 
a blowout jet (iPariat et al.ll2009l) . Following I Yang et alJ (1201 3l) . 
our methods to test the level of magnetic helicity dissipation re¬ 
lies on the comparison of the variation of relative helicity in the 
domain with the fluxes of helicity through the boundaries. 

This lead us to revise the formulation of the time varia¬ 
tion of relative magnetic helicity in a fully bounded volume (cf. 
Sect. 12.3b . As relative magnetic helicity relies on magnetic vec¬ 
tor potential, the question of the gauge is a central problematic 
of any magnetic-helicity related quantity. A general decomposi¬ 
tion of the gauge-invariant time variation of magnetic helicity is 
given in Eq. (|2^ . with no assumption made on the gauges of the 
magnetic held and of the reference potential held. Furthermore, 
we discussed how specific gauges and combination of gauges 
can simplify the formulation of the helicity variation (Sect. lZAl l. 

Following (IValori et al.l 1201 2l) . we computed the variation 
of relative magnetic helicity using different gauges, all of them 
based on the gauge of Eq. (H31 l suggested bv iDe Vo^ (l2000l) . We 
have been able to test the gauge dependance of several terms en¬ 
tering in the decomposition of relative magnetic helicity and its 
time variation. We demonstrated that all the quantities that were 
theoretically gauge-invariant (i/, //j, //pj, d///df) were indeed in¬ 
variant with a very good numerical precision (< 0.3% of relative 
error). 
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Additionally, our analysis showed the effect of using differ¬ 
ent gauges on gauge-dependant quantities. Of particular inter¬ 
est are the results that we obtain on the Fyn and F^n terms, 
Eqs. (126127b . enterin g in the decompo sition of the helicity flux. 
In some studies (e.g. lLiu et al.ll2014bllal) . these terms have been 
used to putatively track the helicity contribution of vertical and 
horizontal plasma flows. However, our computations illustrates 
that these fluxes (and their ratio) vary with the gauges used hence 
precluding any meaningful physic al insight of their interpreta - 
tion in term of helicity. Following iDemoulin & Berg^ (l2003ll . 
we have concluded that only the total helicity flux, Ftot, conveys 
a physical meaning. 

Unlike magnetic helicity, even in ideal MHD, we showed 
that relative magnetic helicity cannot be expressed in general in 
classical conserved form. Only when computing the reference 
potential field with the Coulomb gauge, can the variation of the 
relative magnetic helicity be expressed in ideal MHD as a pure 
surface flux. Generally, relative magnetic helicity is therefore 
not a conserved quantity in a classical sense. 

In the hrst phase of our sim ulation, because of a topologi¬ 
cal constraint (iPariat et al.ll2009h . the system is believed to fol¬ 
low tightly an ideal MHD evolution. De facto, in our practi¬ 
cal numerical case, we observed during that phase, that relative 
magnetic helicity is very well conserved, its variation follow¬ 
ing the time-accumulated flux of helicity with a relative accu¬ 
racy of 0.3%. The relative helicity dissipation that we obtain is 
one order of magnitude smaller than the on e estimated durin g 
the ideal evolution of the simulation tested in lYang et al.l(l2013h . 
The measure of the helicity dissipation can appear as a practical 
way to test the level of ideality in a simulation. During the ideal 
MHD phase, a s expected, the Flux Corrected Transport scheme 
(lDeVorelll991h that is used to produce our test numerical simu¬ 
lation, is effectively able to ensure a quasi-ideal evolution with a 
measurable high degree of precision. 

Furthermore, the term-by-term study of the helicity varia¬ 
tion enables to determine the real, gauge-invariant, dissipation 
of the magnetic helicity of the studied magnetic held d/ Z/dUiss, 
Fa. (l24b . For a solar-like active event, we have conhrmed iTavlorl 
(I I 974 I) hypothesis that magnetic helicity is very well conserved 
even when non-ideal processes are acting. For the specihc event 
that we have studied, even when intense magnetic reconnection 
is present, less than 2.2% of the injected helicity is dissipated. 
While this is one order of magnitude larger than during the ideal 
phase, the dissipation of magnetic helicity is more than 30 times 
smaller than the dissipation of magnetic energy during the same 
per iod. _ 

lYang et ^ (1201 3h and this study are paving the way for fu¬ 
ture, more complete and more extensive tests of the Taylor’s 
conjecture. In parallel to the exploration of the properties of 
helicity, our study offers als o more numericall y-oriented appli¬ 
cations. In a previous work (IValori et al.ll2013h we introduced a 
diagnostic for numerical discretization of magnetic held dataset 
that present hnite error of non-solenoidality which are impact¬ 
ing the estimation of their magnetic energy. Further studies on 
the effect of time and spatial resolution, on a wider range of pro¬ 
cesses and dynamical MHD evolution are now needed. If pre¬ 
cise quantitative bounds are placed on the level of helicity dis¬ 
sipation, magnetic helicity will eventually become useful to test 
numerical MHD codes: the level of helicity dissipation could be 
used as a quantitative criterion of the quality of numerical exper¬ 
iments over the entire simulated evolution. With respect to the 
instantaneous divergence metric of the magnetic held, helicity is 
a complementary, extremely sensitive proxy suitable for testing 
integral conservation properties. The method of the measure of 


the helicity dissipation presented in this manuscript, would open 
up a new way to benchmark numerical codes. 

The strong physical insight that can be gained by studying 
magnetic helicity is also further raised by our present study. 
More that fo rty years after, our numerically precise tests of 
iTavlod (Il974l) conjecture on a solar like events, conhrms that 
magnetic helicity is a quasi-conserved quantity of MHD. The 
application of the conservation of magnetic helicity is full of po¬ 
tential for the study of complex natural and experimental mag¬ 
netised plasma systems. Because of its conservation, magnetic 
helicity may be the raison d’etre o f the existence of coronal 
mass ejection (lRusllll994tlFowlll996h . Magnetic helicity can be 
tracked to characterise and relate the evolution o f coronal active 
region with interplanetary m agnetic clouds (e.g. iMandrini et al.l 
l200llNakwacki et al.ll2011|) . 

Finally, the impact of magnetic helicity conservation on the 
magn etic reconnection m echanism remains to be better under¬ 
stood (iRussell et al.l2015h . While it has been observed that mag- 
neti c helicity can significantly modify the re connection dynam¬ 
ics (iFinton et al.l200l|;lDel Sordo et al.l2010h . how magnetic he¬ 
licity is actually redistributed in the system, at quasi-constant to¬ 
tal value by magnetic reconnection, still needs to be determined. 
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Appendix A: Helicity test without the condition: 

A(ztop) — Ap(ztop) 

In the numerical implementation of both the practical DeVore 
method and the DeVore-Coulomb methods, we impose the con¬ 
dition of Eq. (ESI, of having the same distribution of A and Ap 
at the top boundary, in addition of the DeVore gauge (Eq. l45l l 
which is valid in the whole volume. Because of this condition, 
A and Ap are both different when computed with one method or 
the other. The gauge of A and the gauge of Ap are linked by 
Eq. (l46l l. It induces that Eaap is null at the top boundary, and de 
facto reduces its intensity in the whole domaine 

It is however possible to compute and estimate the helicity 
without enforcing Eq. gai. We can follow the helicity evolu¬ 
tion by mixing the vector potential computed with both method. 
Since gauge invariance of Eq. (fTOl l does not require the use of 
the same gauge for A and Ap, we can use the A computed with 
the practical DeVore method with Ap derived with the DeVore- 
Coulomb methods. Ap is thus still satisfying both the DeVore 
and the Coulomb gauge. As A, which is only satisfying the De¬ 
Vore gauge condition, has been computed independently of Ap, 
there is no boundary surface along which they share any com¬ 
mon distribution. In this appendix we refer to this derivation as 
the "general DeVore-Coulomb" case. 

Eig. IA.il right panel, presents the different terms entering in 
the decompositions Eqs. (II11151) of the relative helicity. As in 
Sect. 15.11 the gauge invariance of in this computation rel¬ 
atively to the other methods is ensured with a high precision 
(<0.3%). As with the others methods, one remarks that //j and 
//pj remains constant while //mix and //p are different in the 
general DeVore-Coulomb case. This further confirms the gauge 
dependance properties of each decomposition. 

When looking at the time integrated helicity fluxes (Eig. IA.il 
right panel), we find again that is following tightly the 

variation of helicity H<v- As for the two other methods, the he¬ 
licity dissipation is also very small and with a precision similar 
to the practical DeVore case (Sect. 1431 ). What significantly dif¬ 
fers from the two other methods is the repartition of the helicity 
flux Ftot between the different terms which composes it. 

Since Ap is respecting the Coulomb condition the term 
d///dfp_var is null to the numerical precision. Since only in¬ 
volves quantities based on the derivation of the potential field, 
Hdv,<i) are equal for both the DeVore-Coulomb and the general 
DeVore-Coulomb cases. On the other hand, as F^n, /^vn are only 
involving A, Flg^ yn and Hgry Bn in the general DeVore-Coulomb 
are equal with their respective estimations in the practical De¬ 
Vore case. 

It is therefore Eaap with concentrate the helicity flux contri¬ 
bution that enable Fioi to be quasi gauge-invariant for the three 
derivations (Eig. lA.ll right panel). While was neg- 

ligeable in both the practical DeVore and the DeVore-Coulomb 
cases, we observe that this term is now a major contributor of the 
helicity fluxes. This is not surprising since Hgry AAp results from 
the existence of large differences between the distribution of A 
and Ap on the boundaries. Both the computations in the prac¬ 
tical DeVore and the DeVore-Coulomb methods were enforcing 
Eq. (l46l l. which induces a very weak value of Hgry AAp- We ob¬ 
serve that dropping the condition (l46l l creates a strong F[g<y^AAp- 

This test further demonstrates that the choice of the gauge 
strongly influences the distribution of the helicity fluxes com¬ 
posing the total helicity flux Etot- Only the total flux Etot is a 
quasi gauge-invariant. None of the terms which are composing 
the helicity flux Etot shall a priori be neglected. Depending on 
the gauge, each term can carry a significant contribution. In a 
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Fig. A.l. Left panel: Relative mag¬ 
netic helicity (Hn/, black dashed line) 
and its decomposition in the general 
DeVore-Coulomb method. The plot¬ 
ted quantities are the same as in Fig¬ 
ures [21 bottom and [7] right panels. 
Right panel: H<v and # evolution 
in the system computed with the general 
DeVore-Coulomb method. The plotted 
quantities are the same as in Fig. [9] 


numerical estimation, it is thus highly advisable to compute all 
the terms which are forming the helicity flux density (Eq. l23l l. 
Explicitly computing each terms allows us to verify that the con¬ 
straints set on the used gauges are effectively enforced numeri¬ 
cally. 
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